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Abstract 

Electron-electron  scattering  conserves  total  momentum  and  does  not  dissi¬ 
pate  momentum  directly  in  a  low-density  system  where  the  Umklapp  process 
is  forbidden.  However,  it  can  still  affect  the  conductance  through  the  energy 
relaxation  of  the  electrons.  We  show  here  that  this  effect  can  be  studied 
with  arbitrary  accuracy  in  a  multi-sublevel  one-dimensional  single  quantum 
wire  system  in  the  presence  of  roughness  and  phonon  scattering  using  a  for¬ 
mally  exact  solution  of  the  Boltzmann  transport  equation.  The  intrasubband 
electron-electron  scattering  is  found  to  yield  no  net  effect  on  the  transport 
of  electrons  in  ID  with  only  one  sublevel  occupied.  For  a  system  with  a 
multi-level  occupation,  however,  we  find  a  significant  effect  of  inter-sublevel 
electron-electron  scattering  on  the  temperature  and  density  dependence  of  the 
resistance  at  low  temperatures. 
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I.  INTRODUCTION 


Carrier  transport  in  low- dimensional  semiconductor  structures  (quantum  wells  and  wires) 
is  of  interest  for  fast  (high  mobility)  electronic  devices.  Many  previous  studies  have  focused 
on  the  low-temperature  conduction  in  semiconductor  quantum  wires. 1-12  The  scope  of  these 
past  investigations  covers  fields  of  low-temperature  clastic  scattering, 1-8  inelastic  scatter¬ 
ing, 10,11  and  magnetic-field  effects 5,6,8-10,12  in  quantum  wires.  In  these  quantum- wire  sys¬ 
tems,  scattering  by  impurities,  interface  roughness  at  low  temperatures  and  by  phonons  at 
high  temperatures  plays  a  dominant  role  for  the  momentum  relaxation  of  drifting  electrons. 
Although  electron-electron  scattering  has  been  known  to  affect  the  carrier  transport  in  semi¬ 
conductors  for  a  long  time, 13-15  it  has  not  received  sufficient  attention  in  the  past,  especially 
in  the  recent  calculations  of  mobilities  in  semiconductor  quantum-wire  systems,  due  to  the 
complexity  of  calculations  involved.  In  this  paper,  we  examine  the  effect  of  electron-electron 
scattering  on  the  resistivity  using  the  Fermi  liquid  model.  As  is  well  known,  electron-electron 
interaction  yields  the  Liittinger  liquid  effect  in  clean  systems.  This  effect  will  not  be  con¬ 
sidered  here. 

In  modulation  n-doped  quantum-wire  structures,  there  are  usually  many  more  conduc¬ 
tion  electrons  than  ionized  impurities  inside  the  channel.  As  a  result,  an  electron  can  be 
scattered  much  more  frequently  by  other  electrons  than  by  ionized  impurities,  leading  to 
an  enhanced  effect  of  electron-electron  scattering.  In  addition,  acoustic  phonon  scattering 
is  suppressed  at  low  temperatures. 10,11  Therefore,  electron-electron  scattering  cannot  be  ig¬ 
nored  in  the  calculation  of  the  conductance  in  modulation-doped  quantum  wires.  However, 
electron-electron  scattering  alone  does  not  produce  momentum  dissipation  directly  because 
of  two-particle  momentum  conservation. 16  The  major  role  of  electron-electron  scattering  is 
the  energy  relaxation  where  carriers  undergo  rapid  inelastic  transitions  among  all  the  states 
accessible  by  electron-electron  scattering.  As  a  result,  the  momentum  relaxation  rate  be¬ 
comes  averaged  over  many  states,  becoming  less  dependent  on  the  energy,  and  is  driven 
from  the  original  exact  relaxation  rate  for  the  rest  of  the  elastic  and  inelastic  scattering 
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mechanisms  without  electron-electron  scattering.  This  effect  increases  the  resistivity.  A 


quantitative  argument  based  on  the  variational  principle  was  given  earlier  for  this  effect.15 
In  this  paper,  we  quantitatively  study  the  effect  of  the  competition  of  elastic,  inelastic  and 
electron-electron  scattering  as  functions  of  the  temperature  T  and  the  electron  density  nn> 


The  intrasubband  electron-electron  scattering  is  found  to  yield  no  net  effect  on  the  trans¬ 
port  of  electrons  in  the  ID  limit  where  only  one  sublevel  is  occupied.  This  is  in  sharp 
contrast  to  a  finite  intrasubband  electron-electron  scattering  in  the  2D  limit. 15  In  multi¬ 
sublevel  quantum  wires,  intersubband  electron-electron  scattering  is  found  to  be  significant 
at  low  temperatures  in  momentum-relaxation  processes  as  will  be  demonstrated  later  in  our 
numerical  calculations.  It  is  prohibitively  difficult  to  study  the  effect  of  electron-electron 
scattering  in  two  and  three  dimensions  accurately.  In  a  one  dimensional  system,  however, 
we  can  develop  an  accurate  solution  to  the  Boltzmann  transport  equation  utilizing  the  fact 
that  there  are  only  a  discrete  number  of  points  in  k  space  for  a  given  energy  unlike  in  higher 
dimensions. 

We  establish  a  matrix-equation  approach  for  numerical  calculation.  This  approach  has 
the  advantage  of  performing  numerical  calculation  with  any  desired  accuracy.  As  a  result, 
the  Boltzmann  transport  equation  (described  in  the  next  section)  can  be  solved  exactly 
formally  for  quantum  wires  yielding  accurate  numerical  results.  We  apply  this  approach 
to  study  the  temperature-  and  density-dependent  resistance  of  a  single  quantum  wire.  For 
elastic  scattering,  we  consider  only  surface-roughness  scattering  and  neglect  scattering  from 
ionized  impurities  by  assuming  that  dopants  sit  far  away  from  the  quantum  wire  in  the 
growth  (z)  direction  in  modulation-doped  systems.  We  study  only  high-quality  quantum 
wires  where  the  localization  length  is  expected  to  be  longer  than  the  sample  length. 

The  outline  of  this  paper  is  as  follows.  In  Sec.  II  we  derive  an  expression  for  the  con¬ 
ductance  in  the  presence  of  interface-roughness,  phonon  and  electron-electron  scattering  in 
a  multi-sublevel  single  quantum  wire.  In  Sec.  Ill  we  construct  a  matrix-equation  approach 
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for  solving  the  Boltzmann  transport  equation  numerically  with  any  desired  accuracy.  In 
Sec.  IV  effects  of  intersubband  electron-electron  scattering  on  the  temperature  and  density 
dependence  of  the  resistance  are  computed  numerically,  along  with  detailed  discussions  and 
explanations  of  the  numerical  results.  A  brief  conclusion  is  given  in  Sec.  V. 

II.  SCATTERING  EFFECTS  ON  MULTI-SUBLEVEL  TRANSPORT 

In  this  section,  we  briefly  review  the  key  results  in  our  previous  studies8,10  and  repeat 
some  of  the  equations  for  completeness  of  the  theory  presented  in  this  paper.  By  using 
notation  nearly  identical  to  that  already  in  the  literatures, 8,10  the  Boltzmann  transport 
equation  for  electrons  along  the  wire  (y)  direction  is  given  by 

2tt 

vj  +  —  ^2  \Ij',j\  idj'  —  9j )  8  {£j  ~~  £j')  +  Pj  +  Qj  =  0  ,  (1) 

i' 

where  Pj  stands  for  the  electron-phonon  scattering  and  Qj  represents  electron-electron  scat¬ 
tering.  These  quantities  will  be  given  below.  In  Eq.  (1),  £j  and  Vj  =  U~ld,£j/dk  are  the 
kinetic  energy  and  the  group  velocity  of  the  electrons  in  the  electronic  state  j  =  [n,  k}. 
Here,  n  —  0,  1,  •  •  •  represents  the  sublevel  index  due  to  the  size  quantization  of  the  wire  in 
the  x  direction  and  k  is  the  wave  number  of  free  electrons  along  the  y  direction.  Moreover, 
Ijtj  is  the  elastic  scattering  matrix8  clue  to  interface  roughness.  The  quantity  gj  is  pro¬ 
portional  to  the  nonequilibrium  deviation  from  the  equilibrium  Fermi  distribution  function 
fj  ^  =  fo(£j)  after  a  linear  expansion  of  the  distribution  function  to  the  leading  order  in  the 
external  electric  held  E :  fj  —  fj  ^  +  g j[—d fo{£ j)de j]eE . 

The  quantity  Pj  in  Eq.  (1),  representing  the  contribution  from  the  inelastic  electron- 
phonon  interaction,  is  given  by10,1' 

Pj  =  Y  Y,  I vj'9j\2  (/j'T)  +  nsq)  ( 9j '  -  9j )  8  (£j  -  £f  ±  husf)  5kfy  k±qy  ,  (2) 

j',s,q,± 

where  V^j  is  the  screened  electron-phonon  interaction 11  dehned  by  <  ri'k’  nk  >= 

VPj  5k \  k+qy  and  the  sign  —  (or  +)  corresponds  to  the  one-phonon  emission  (or  absorption) 
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process.  In  Eq.  (2),  /j,+)  =  1  —  /■,  J  and  ns^  is  the  Boson  distribution  function  for  phonons 
of  mode  s ,  having  wave  vector  q  =  (qx,  qy,  qz)  and  equilibrium  phonon  energy  hcos 

To  consider  the  effect  of  electron-electron  scattering  on  transport  we  assume  that  only 
the  quantum-well  ground  state  <j>o(z)  with  energy  level  Eiz  is  relevant  in  the  z  direction  due 
to  an  extra-thin  thickness  of  the  wire  in  this  direction  and  low  temperatures.  The  wave 
function  along  the  wire  direction  is  just  a  plane  wave.  The  wave  function  of  the  nth  sublevel 
in  the  x  direction  is  denoted  by  ij ’n(x ).  We  further  assume  that  the  confining  potential  in 
the  x  direction  is  a  parabolic  confinement.  This  gives  rise  to  a  harmonic  wave  function 
—  (a/v^h^n!)1/2  exp(— a2x2 /2)Hn(ax)  and  a  quantized  energy  level  (n  +  1/2)HQq, 
where  Hn(x)  is  the  nth-order  Hermite  polynomial,  a  =  ^m*£lo/h,  m*  is  the  effective  mass 
of  electrons,  and  M70  is  the  energy-level  separation.  An  electron  in  an  initial  state  j  =  (n,  k) 
will  be  scattered  into  a  final  state  j'  =  (n',  k')  by  simultaneously  scattering  another  electron 
from  an  initial  state  j\  =  (n i,  k\)  into  a  hnal  state  j[  =  {n\ ,  k[)  through  the  Coulomb 
interaction  between  two  electrons,  which  is  given  by 

Klh  =  — “72  /  rf3'r  ^o(-)  C'O^nO)  exp  [i  (k  -  k')  y] 

€q erLz  J 


x  /  d3ri  il>*fi(xi)ij>ni(xi) 


exp  [i  (A?!  -  k[)  yi] 


|  r  —  rq  | 

where  er  is  the  dielectric  constant  of  the  host  material  and  L  is  the  wire  length.  This  can 
be  explicitly  calculated  as 

02  r  p 


3'n  e0evL 


$ qy ,  ki—k'^k+ki,  k'+k^  dx  dz  1p^,(x)tpn(x) 


x  /  dXl  I dZl  ^ Zl )  /  dv  ^r==f  > 

and  is  further  simplified  to 


/v'V:  =  5, 


,  (Jfl) 

\e0erL) 
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x  J  dx  J  dz  J  dx\  J  dzi  K0(\k'  -  k\TZ)  <j>l(z)il)*l(x)i/)n(x)  <f>0(z1)2^(x1)i/}  ni(ah)  ,  (3) 

where  K0(x)  is  the  modified  Bessel  function  of  the  third  kind,  and  7 Z  = 

\J (zi  —  z)2  +  (aq  —  x)2 .  By  denoting  the  two-electron  wave  function  as 

.T,  /  •  •  \  ,  /  ^exp  (iky)  exp  (ikm) 

J l)  =  ^n{x) - ~j= -  ^ni(Zl) - ~j= -  , 

the  matrix  element  of  the  Coulomb  interaction  Vc  for  the  normalized  spin  singlet  (upper 
sign)  and  the  spin  triplet  (lower  sign)  states  is  given  by 

l  <  ±  Wi,f)  IKI  ±  > 


=<  \  vc\  vh(i, Jd)  >  ±  <  |K|  *(j, iO  >  • 

Here,  the  first  term  corresponds  to  Eq.  (3)  and  the  second  (he.,  the  exchange)  term  is  found 
by  interchanging  j'  j[.  The  interaction  averaged  over  the  spin  configuration  leads  to 


K =  KJ 

3,31  3,3 1 


f’fi  1  t sii’i’ 


K  1 

2  j’jl 


(4) 


The  electron-electron  scattering  adds  a  term  Qj  on  the  left  hand  side  of  Eq.  (I),18,19  which 
is  given  by 


Qj  = 


47 r 

hk^T 


df, 


(-) 


dej 


-i 


E 

j'ju'i 


701 


f  (-)  f  H  f  (+)  f  (+) 

Jj  Jjl  J  j'  J  j[ 


x  (&'  +  9j[  -  9j  -  9ji)  ^  {ej  + 


£j  i  Ah 


(5) 


Here,  T  is  the  temperature  of  the  thermal-equilibrium  system,  and  the  spin  degeneracy  has 
been  taken  into  account.  The  dynamical  screening  to  the  electron-electron  scattering  in 
Eq.  (5)  can  be  included  through  the  following  substitutions 


/  Kr’l[ 

\  3,3 1  ’  3,3 1  / 


3',3x 


K  : 

3,3 1 


K: 


3,31 


efjA[\k>  ~  k\i  ~  £j)/h ]  ’  -  kl  (£j[  -  £j)/n\ 

where  ePPA ( | qy | ,  oj)  can  be  found  from  diagonal  elements  of  a  dielectric-function  matrix 
under  a  random-phase  approximation  (RPA).20  After  the  Fourier  transform  of  K0(x)  with 

jf  jf 

respect  to  rj,  the  expression  for  ^-dependent  Kj^1  in  Eq.  (3)  is  found  to  be 
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^  ^n>—  n<,even  dn>  _n<  eve  n  2(5n>  _n< ;  0(jd  n<,odd)  1 

where  qy  =  k'  —  k,  L^\x)  is  the  nth-order  associated  Laguerre  polynomial,  n<  =  min(n,  nr) , 
n>  =  max(n,  n'),  nf  =  min(ni,  n\ ) ,  nf  =  max(ni,  n^),  P  =  rj2 /2a2,  and  the  form  factor 

is 


9i/)  =  f  dz  J dzi  ^o(-)  exp(-\/h2  +  ^k  -  zi\)  Mz i)2  • 


The  conductance  is  found  to  be8 

2  _  /*  +  00  /Q 

G  =  JJ  J2  vj9j  Jo  de  5  (, £j  -  e)  -^/(_)(£) 

2p2  /*+oo  f)  /  \ 

=  hiL  dE  [~Tef  (£)J  5 (6) 

where  the  symbol  (8)  represents  the  matrix  inner  product,  the  /e-summation  for  the  running 
variable  j  is  converted  into  the  ^-integration  in  the  second  equality,  and  S,  g  are  column 
vectors  with  elements  sv  and  gu  defined  below.  The  index  v  describes  a  group  of  discrete 
equi-energy  points  given  by  the  intersection  of  the  horizontal  line  at  an  energy  £  with  a 
series  of  energy-dispersion  curves,  namely  the  solution  of  £  =  e3  for  each  j.  The  element 
su  =  vv/  \vu\  of  the  column  vector  S  in  Eq.  (6)  is  the  sign  of  the  slope  of  the  tangent  to  the 
dispersion  curves  at  these  points.  The  quantity  gu  is  related  to  gu  at  the  equi-energy  point 
v.  Enumerating  the  indices  of  these  points  by 

v  = —N,  —N  +  l,  •••,  —1,  1,  •••,  N  —  l,  N  (7) 
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for  symmetric  dispersion  curves,  introducing  gu  =  gu+gN,  and  using  the  fact  that  g_v  =  — 
we  find8 

9v  =  %—  2  9-n,  9n  =  0.  (8) 


The  branch  index  n  and  the  wave  number  k  of  the  equi-energy  point  v  depends  on 
the  energy  e.  Hereafter,  we  put  this  e-dependence  specifically  on  all  quantities,  e.g.,  gv(e). 
Following  the  previous  treatment,8  we  find 

U{e)  ®  g{e)  +  V_1(e)  ®  P{e)  +  V_1(e)  (8)  Q(e)  =  -S(s)  .  (9) 


Here  g,  S,  P  and  Q  are  column  vectors  with  the  index  v  excluding  the  end  point  u  =  N 
for  each  energy  e  in  Eq.  (7),  V(e)  is  a  diagonal  matrix  with  absolute  velocity  |u|  being  its 
element  in  u,  and  U  is  the  2 N  —  1  by  2N  —  1  square  matrix  (taking  away  the  last  row  and 
column)  8  with  its  off-diagonal  elements  (j'  ^  j )  by 


uj',j  ~ 


aA2V^5b2 

h2\yj'vj\ 


exp 


-(k  -  k')2K2 


\MziW 


r+o o 

x  /  dg  exp(-r)2/2)Ln(rj2/2)Ln/(r]2/2) 

Jo 

and  its  diagonal  elements  by 

uj,j  =  ~  uhi'  ■ 
iVi 

Here,  Ln(x)  =  L^(x),  Z\  is  the  position  of  the  AlrGai^xAs/GaAs-interface  of  a  symmetric 
quantum  well  in  the  growth  direction,  Vo  is  the  conduction  band  offset  at  the  interface,  6b 
is  the  average  layer  fluctuation,  and  A  is  the  correlation  length  of  roughness  in  the  x,  y 
directions.  The  elimination  of  the  TVth-point  in  the  process  of  reducing  Eq.  (1)  into  Eq.  (9) 
is  necessary  clue  to  the  fact  that  the  matrix  U  is  singular  in  general  and  has  no  inverse  if 
this  point  is  included.  8 


In  order  to  carry  out  the  ('/-summation  in  Eq.  (2)  for  phonon  scattering,  we  introduce  a 
cylindrical  coordinate  qx  =  q±  cos  (j),  qz  =  q±  sin  4>,  q  =  \Jq\  +  and  transform 

Y,  d  Qy >  k'-k  =  ^27t)2  J  f**  ^1y’  k'~k  =  ^27r)2  j  J  d4>  Sqy,  k'-k  , 

where  A  is  the  cross-sectional  area  of  the  wire.  We  now  employ  the  Debye  approximation 
at  low  temperatures  for  acoustic  phonons  and  use  Hus$  =  hcsq,  where  cs  is  the  sound 
velocity.  Carrying  out  the  ^-integration  through  the  energy  delta  function  and  converting 
the  //-summation  into  the  energy-integration,  we  find 


Pu(e) 


/(e)  |  (2vr)Vnf 


Y  d£ 


,  <?©  (±e'  =F  e) 


V&(e,e') 

\vj'{e')vv{e)\ 


x  (jf  +  ns2)  [gAe')  -gv{e)\  , 


where  fl  =  AL  is  the  sample  volume,  ©(e)  is  the  unit  step  function, 


hcsq  =  ev  -  e',  ,  \qy\  =  \ k'  -  k |  ,  and  q±  =  J q2  -  q2 


In  Eq.  (10),  eu  indicates  the  point  v  at  the  energy  e,  and  the  s'- integration  excludes  the 
region  where  q±  becomes  imaginary.  The  electron-phonon  interaction  matrix  elements  in 


Eq.  (10)  are  calculated  as 


l^n'fnl  —  I  RPA/'I  J\  0  12  ^n'n (qx) Az (?z) 

I erin  {IQyU^jl 


with  u>  =  (sn/jk,  -  £n,k)/h,  where 


A«'n(fe)  =  -^iTexp APAL^iP,)}2  , 
n>- 

Pi  =  ql/2a2,  m  =  n>  —  n<,  n>  (n<)  are  the  larger  (lesser)  of  n  and  n',  and 


A z(qz)  =  J  dz  exp {iqzz)4>l(z) 


Moreover,  for  the  longitudinal  (s  —  l)  and  transverse  (s  =  t)  modes  of  acoustic  phonons,  we 
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1^1=  — 

lq  2QpoCl 


D2  +  (e/114) 


2  Mo) 


v.i  = 

q  2Op0  ct 


kq  ,_L 


(1 eh 
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where  po,  D ,  and  hu  are  the  ion-mass  density,  deformation-potential  coefficient,  and  piezo¬ 
electric  constant 10  and  the  form  factors  are 


M9)  = 


Mi)  = 


2[q2(qlq2y  +  q2yq2z  +  q2q2x)  -  9 q2q2qz. 


By  defining 

h2f 

£j(k  =  0)  +  eh{k  1  =  0)  -  £f(k'  =  0)  -  £fi(k[  =  0)  =  —  , 

TYl 

the  energy-delta  function  in  Eq.  (5)  with  a  parabolic  energy  dispersion  £j  =  =  Eiz  + 

(n  +  l/2)hklo  +  h2k2/2m*  can  be  rewritten  as 

5  (ej  +  £h  -  £f  -  £jj)  =  ^  5  P'  ~  k)  (k'  -  ki)  -  f]  . 

Here,  C,  is  related  to  the  energy-level  separation  M20.  When  both  electrons  undergo  intra¬ 
subband  scattering,  only  the  exchange  scattering  with  k!  =  k\  and  k[  —  k  is  possible  in 
view  of  £  =  0.  In  this  case,  however,  Qj  =  0  because  the  sum  of  all  g  factors  in  Eq.  (5) 
vanishes.  Therefore,  the  intrasubband  electron-electron  scattering  cancels  itself  completely 
and  has  no  net  effect  on  the  conductance  for  a  single-sublevel  occupation  in  the  ID  limit. 
In  the  following,  we  consider  intersubband  electron-electron  scattering  that  involves  only 
two  occupied  sublevels  at  low  densities  for  simplicity  and  write  HQ0  =  h2k\/2m*.  Possible 
values  of  £  are  given  by  £  =  5nk\/ 2,  where  Sn  =  ±1,  ±2,  depending  on  whether  only  one 
electron  undergoes  intersubband  scattering  or  both  electrons  jump  from  the  same  sublevel 
to  the  other. 
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For  electron-electron  scattering,  we  now  carry  out  momentum  summations  in  Eq.  (5)  for 
the  last  two  terms  oc  g3  +  gn  which  are  denoted  by  Qj  \  For  this  purpose,  we  write  the 
energy  delta  function  as 

stc 

8  (ej  +  £jx  —  £ji  —  £j'^j  =  -j-j-  8  ([k1  —  k_(ki,  k,  8n )]  [k'  —  k+(ki,  k,  <5n)]) 


m*0  [D{ki  —  k,  8n )] 
h2 yj D{k\  —  k ,  8n) 


{<5  [k'  —  &_(&!,  k,  5n)\  +  8[k'  —  k+(ki,  k,  5n)]}  , 


where 


k±(kl,k,8n)  =  - 


k1  +  k  ±  \J D(ki  —  k,  5n) 


,  D(k\  —  k,  8n)  =  (£q  —  k)2  +  28nk\  ,  (12) 


8n  —  n  +  n\  —  n'  —  n[  ^  0.  From  these,  we  hnd 

777*  L2 


QP  = 
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3  (  |  ^ 

fj  ni,n',nj,s=di 


dk\ 


isi'j 1 12 
3,3 1 


xtfW  te  +  fe)  9P(<:i“<:'fa)1 


\jD(k\  —  k,  8n) 


8k j,  ks(ki,k,Sn)  8k',  k—s(ki,k,Sn)  ,  8 Tl  7^  0  .  (13) 


The  first  two  terms  oc  gy  +  g}^  in  the  momentum  summations  of  Eq.  (5)  are  denoted  by 
Q(j+)  ■  Using 

8  ( £j  +  £j1  —  Eji  —  Ey^j  =  &  ([ki  —  k-(kr,  k[,  —  8n)]  [k\  —  k+(k',  k[,  —  8n)]) 


777*0  [D{k[  —  k',  —  8n)] 


{5  [ki  —  k_(kr,  k[,  —8n)]  +  8  [ki  —  k+(k',  k[,  —  <5 77) ] }  , 


h2yjD(k[  -  k\  —8n) 

we  hnd  k\  =  k±(k',  k\ ,  —8n)  and  k  =  kT(k',  k\ ,  — 8n ).  The  expression  for  k  yields 


iJ  -  u _ _ 

Ki~  K  oth, 


8nk\  ,  , .  ,  8nk‘i 


2  (k’  -  k) 


k'  =  k- 
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Based  on  these,  we  arrive  at 
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We  now  interchange  variables  j'  < — »  j\  for  the  ^'-integration  and  variables  < — >  j\  for 

the  /^-integration  in  Eq.  (15),  and  obtain 
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JJi 


'J1  ^  ^  \/D{ki  -  fc',  -<Jn2) 
where  Sn j  =  n  —  ri\  +  n'  —  n'j  7^  0  and  5n2  =  n  —  rq  —  n'  +  n[  ^  0. 


III.  MATRIX  FORMALISM  FOR  THE  CONDUCTANCE 


In  this  section,  we  present  details  of  the  matrix  formalism  for  solving  the  Boltzmann 
transport  equation  numerically  with  any  desired  accuracy.  For  the  phonon  scattering  P 
term,  let  us  first  define  the  scattering-out  (second)  term  on  the  right-hand  side  of  Eq.  (10) 
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de 


,  <?9  (±e;  =F  g) 
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d(j> 


VAMe')  ,  (T)  , 

iv(£)iy(£')l  ^  1  7) 
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n' s  ±  u  ,VL'S  °  I 

along  with  conditions  given  in  Eq.  (11).  We  now  replace  the  ^'-integration  in  the  scattering- 
in  (first)  term  in  Eq.  (10)  by  a  summation  over  energy  points  uniformly  spaced  with  a 
sufficiently  small  energy  interval  Se.  We  further  introduce  a  new  index  t  =  {£,  m},  where  £ 
counts  points  from  the  left  to  the  right  for  a  given  energy  and  m  indicates  energy  branches 
in  the  above  subdivision.  Defining 


Wtiif  = 


f 16 e  ^  qQ  (±£7  =F  e) 

~2  E 


(2ir)2h2 


S,d= 


he. 


/  #  44  (f™  +  n«)  ’ 


(18) 
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we  can  recast  Eq.  (10)  into  a  matrix  form 


A  +  w}o)  gt  =  J2  Wtj  gt>  (19) 

with  Hyoj  =  J2t'  This  equation  can  be  rewritten  in  view  of  Eq.  (8)  as 

R  =  ?  *•'  *■  (20) 

where  for  t  =  {£,  m}  and  t'  =  {£' ,  m'}  we  have 

Wt>?  =  ( Wt,t ’  -  St,?  Wt{o))  (l  -  -  .  (21) 

It  is  understood  that  the  index  t  (t')  excludes  the  right-most  point  for  a  given  energy  index 
m  (m').  Here,  £’  =  1  signifies  the  left-most  point  for  any  energy  index  wl . 


For  the  intersubband  electron-electron  scattering  Q  term,  we  now  transform  the  k\- 
integration  in  Eqs.  (13)  and  (16)  into  an  energy  integration  using  dk\  =  d£j1/H  \vJl  |  and 
change  the  r, -integration  into  a  discrete  summation  by  chopping  the  integration  into  many 
sufficiently  small  uniform  intervals  of  width  Se  as  before.  Equation  (13)  leads  to 
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The  quantity  Q E  in  Eq.  (16)  can  also  be  written  in  a  matrix  form 
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(24) 


where 


Inserting  Eqs.  (20)  and  (27)  into  Eq.  (9),  we  find 

U  +  W  +  z]  (8)  g  =  -S  ,  (29) 

and  the  solution  of  Eq.  (29)  is 

g  =  -  [U  +  W  +  Z]'1  ®  S  .  (30) 

Here,  indices  t  =  {£,  m}  and  t'  =  {£',  rn'}  exclude  the  right-most  points  for  any  given  energy 
indices  m,  m! .  The  index  m  is  the  k  value  at  the  intersection  point  t.  Note  that  U  is  block- 
diagonal,  namely  Ut^  oc  Sm^m,  while  \vt^\  oc  Stjt>  is  diagonal.  When  both  the  roughness  and 
phonon  scattering  are  absent  from  Eq.  (30),  i.e.  U  =  W  =  0,  we  find  that  the  matrix  Z  itself 
is  singular.  This  confirms  that  the  electron- electron  scattering  alone  does  not  contribute  to 
the  resistance  in  the  transport  due  to  the  fact  that  it  conserves  total  momentum  and  does 
not  dissipate  momentum.  The  conductance  is  found  from  the  discrete  version  of  Eq.  (6): 
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IV.  NUMERICAL  RESULTS  AND  DISCUSSION 


In  this  section,  we  present  the  numerical  results  and  discussions.  In  our  numerical  calcu¬ 
lations,  we  have  chosen  the  following  parameters.  For  the  quantum  well  in  the  z  direction, 
we  employ  m*  =  0.067  m0  in  the  well  (m0  being  the  free  electron  mass),  tub  =  0.073  m0 
in  the  barrier,  Vo  =  280  rneV  for  the  barrier  height,  Lw  =  50  A  for  the  well  width,  and 
er  =  12  for  the  average  dielectric  constant.  These  parameters  yield  Elz  =  84  rneV  for  the 
ground-state  energy  and  E2z  —  Elz  =  188  rneV  for  the  energy  separation  between  the  first 
excited  level  and  the  ground  state.  For  the  parabolic  confinement  in  the  x  direction,  we 
use  hflo  =  2.7  rneV  for  the  energy-level  separation.  For  the  interface  roughness,  we  as¬ 
sume  5b  =  5  A  for  the  average  layer  fluctuation  in  the  z  direction  and  A  =  10  A  for  the 
Gaussian  correlation  length  in  the  x,  y  directions.  For  the  acoustic  phonons  in  the  bulk, 
the  following  parameters  are  employed:  Cj  =  5.14  x  105  cm/s  for  the  longitudinal  mode, 
Ct  =  3.04  x  105  cm/s  for  the  transverse  mode,  D  =  —9.3  eV  for  the  deformation  potential, 
hu  =  1.2  x  10'  V ) cm  for  the  piezoelectric  field,  and  p0  =  5.3  g/cm 3  for  the  ion  mass  density. 
The  wire  length  in  the  y  direction  is  assumed  to  be  L  =  1  ym  except  for  the  cases  indicated. 

For  the  ID  density  riio  =  106  cm-1  assumed  for  the  numerical  study,  the  relative  chemical 
potential  y  =  y  —  Elz  —  M20/2  shows  that  M20  <  y  <  2hfl0  for  the  temperature  range 
0.5  <  T  <  8  K  considered,  indicating  occupation  of  the  two  lowest  sublevels  n  =  0,  1. 
When  T  tends  to  zero,  y  approaches  the  Fermi  levels  E1F  =  4.98  rneV  and  E2F  =  2.28  meV, 
indicating  that  the  electrons  are  highly  degenerate. 

In  Fig.  1,  we  present  the  scaled  resistances  R/Ro  as  a  function  of  the  temperature  T  in 
the  range  4  <  T  <8  K.  Here,  Rq  is  the  resistance  R\  due  to  surface  roughness  scattering 
at  T  —  0.5  K .  When  only  the  roughness  scattering  is  present,  R\  (dash-dot-dotted  curve) 
increases  with  T  due  to  gradual  population  ( y  <  2 M20)  of  the  bottom  of  the  sublevel  n  =  2, 
where  the  density-of-states  (proportional  to  the  inverse  of  the  electron  group  velocity)  is 
very  large,  yielding  enhanced  roughness  scattering.  When  T  is  very  high,  R\  becomes  linear 
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in  T.  When  phonon  scattering  is  added,  the  total  resistance  R2  (dashed  curve)  increases 
with  T  faster  compared  to  R\ ,  causing  an  appreciable  enhancement  of  the  resistance  at 
T  —  8  K.  When  electron-electron  scattering  is  further  included,  i?3  (solid  curve)  displays 
a  much  stronger  increase  with  T  compared  to  R2,  showing  a  significant  enhancement  of 
the  resistance  at  T  —  8  K .  On  the  other  hand,  the  differences  among  R,\ ,  R2,  and  i?3 
become  negligible  at  low  temperatures  T  <  1  K .  The  inset  of  the  figure  shows  R3/R0  as  a 
function  of  T  in  the  whole  range  of  T  considered.  The  overall  feature  of  R3/R0  reflects  the 
T-dependence  of  R1/R0  since  the  system  is  dominated  by  roughness  scattering. 

Electron-electron  scattering  itself  does  not  directly  contribute  to  the  resistance  in  the 
transport  of  electrons  because  it  conserves  the  two-particle  momentum  through  the  scat¬ 
tering  process.  However,  it  contributes  to  the  momentum  dissipation  indirectly  by  redis¬ 
tributing  the  electron  energy.  This  point  can  be  made  clear  by  using  Kohler’s  variational 
principle. 15,18  According  to  this  theorem,  the  resistance  is  minimum  and  exact  when  the 
entropy  production  rate  due  to  the  external  field  is  minimum.  To  make  the  argument  sim¬ 
ple,  we  suppose  for  a  moment  that  the  electron-electron  interaction  is  infinitely  large  in 
comparison  with  other  elastic  and  inelastic  interactions.  In  this  case,  the  exact  resistivity  is 
obtained  when  the  transport  relaxation  time  is  independent  of  the  energy  in  the  effective- 
mass  model.  The  exact  resistivity  is  proportional  to  the  entropy  production  rate  under  the 
DC  field  which  vanishes  for  this  solution.18  What  happens  is  that  the  electrons  are  swept 
through  all  possible  accessible  states  rapidly  by  electron-electron  scattering.  As  a  result, 
the  relaxation  rate  is  averaged  over  these  states  and  becomes  independent  of  the  energy. 
While  this  energy-independent  relaxation  rate  minimizes  the  entropy  production  rate  from 
electron-electron  scattering  and  the  resistance,  this  solution  does  not  minimize  the  entropy 
production  rate  from  other  elastic  and  inelastic  scattering  mechanisms  and  raises  the  resis¬ 
tivity  for  these  mechanisms  beyond  their  original  exact  resistivity  obtained  in  the  absence 
of  electron-electron  scattering.  In  the  presence  of  realistic  (he.,  finite)  electron-electron 
scattering,  there  is  a  compromise  between  the  electron-electron  scattering  and  rest  of  the 
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scattering  mechanisms  in  minimizing  the  entropy  production  rate,  enhancing  the  net  re¬ 
sistivity  beyond  the  resistivity  value  without  electron-electron  scattering.  Practically,  the 
roughness  scattering  often  dominates  the  phonon  and  electron-electron  scattering  at  low 
temperatures,  as  seen  in  Fig.  1.  However,  we  can  theoretically  highlight  the  importance  of 
the  buried  electron-electron  scattering  either  by  deducting  the  contribution  of  roughness 
scattering  or  even  by  assuming  the  absence  of  roughness  scattering  in  our  system. 

Figure  2  shows  the  effect  of  electron-electron  scattering  on  the  resistivity  in  the  presence 
of  strong  roughness  scattering  in  (a)  and  without  roughness  scattering  in  (b).  Case  (a) 
corresponds  to  the  situation  studied  in  Fig.  1  with  the  same  definitions  for  Rly  R2,  R3 
but  not  for  i?o-  The  contribution  from  the  phonon  or  electron-electron  scattering  can  be 
separated  as  clearly  seen  in  (a)  which  presents  the  scaled  resistance  differences  A R/Rq 
and  in  (b)  which  shows  the  scaled  resistance  R/Ro  as  a  function  of  the  temperature.  In 
(a),  the  contribution  from  phonon  scattering  is  shown  as  R-2  —  R\  (dashed  curve),  while  the 
contribution  from  electron-electron  scattering  is  shown  as  R3  —  R2  (solid  curve).  The  phonon 
scattering  part  R,2  —  R,\  displays  a  rapid  rise  above  T  >  T0  =  2cs^2m*(p,  —  M2o)Ab  ~  2.7  K. 
The  energy  /crTq  equals  the  energy  of  the  phonon  with  wave  number  qy  =  2k±p  corresponding 
to  the  momentum  transfer  between  the  Fermi  points  in  the  sublevel  n  —  1.  For  T  T0, 
R-2  —  R\  is  found  to  exhibit  a  linear  behavior  in  T/Tq.  The  electron-electron  scattering 
contribution  R3  —  R2  vanishes  at  T  =  0  K  due  to  the  restricted  phase  space  available 
for  scattering.  However,  R3  —  R-2  increases  rapidly  with  T,  roughly  oc  T2  below  2  K , 
a  well-known  behavior  from  the  Umklapp  scattering  process.18  In  the  high-temperature 
nondegenerate  regime,  we  can  show  analytically  that  ( R3  —  R2)/R\  approaches  a  constant 
(4  —  7r) / 7r  for  the  roughness-scattering  dominated  system  in  the  presence  of  strong  electron- 
electron  scattering  following  the  method  given  earlier  by  one  of  the  authors.15  In  this  high 
temperature  regime,  the  increase  of  R3  —  R-2  slows  down  to  a  quasilinear  dependence  of  R\  as 
seen  from  (a).  In  contrast,  we  do  not  see  the  slow  down  in  the  increase  of  electron-electron 
scattering  contribution  with  T  at  high  temperatures  in  (b)  in  the  absence  of  roughness 
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scattering.  In  this  case,  we  can  show  analytically  that  the  ratio  of  the  upper  solid  curve  and 
the  lower  dashed  curve  approaches  a  constant  4/7T  in  the  high-temperature  nondegenerate 
regime  for  strong  electron-electron  interaction.15  Interestingly,  the  numerical  ratio  of  the 
resistance  between  the  solid  and  dashed  curves  in  (b)  approaches  roughly  4/7T  ~  1.3  at  high 
temperatures. 

In  the  following  we  study  the  density  dependence  of  the  resistivity  at  T  —  4  K.  The 
quantity  ft  increases  with  uid  from  one  range  0  <  ft  <  to  another  range  <  ft  <  2Mlo, 
passing  through  the  occupation  of  the  bottom  of  sublevel  n  —  1  in  the  x  direction  at 
n i d  =  0.5  x  106  cm _1.  In  Fig.  3,  we  show  the  calculated  scaled  resistances  R/ Rq  as  functions 
of  the  density  n \ d  within  a  small  range  0.8  x  106  <  n \ d  <  1.0  x  106  cm ~1.  Here,  Rq  is  the 
resistance  Ri  due  to  surface  roughness  scattering  at  nix>  =  106  cm-1.  When  only  roughness 
scattering  is  present,  the  resistance  R\  (dash-dot-dotted  curve)  decreases  with  n \ d  as  I/tiid- 
With  phonon  scattering  added,  we  see  very  small  changes  in  R2  (dashed  curve).  When 
electron-electron  scattering  is  further  taken  into  account,  we  find  relatively  large  changes  in 
f?3  (solid  curve)  compared  to  that  in  R2.  In  addition,  the  differences  among  Ri,  R2,  and  R% 
as  functions  of  Uid  are  smaller  compared  to  those  shown  as  functions  of  T  in  Fig.  1.  The  inset 
of  the  figure  shows  R3/R0  as  a  function  of  n id  in  the  whole  range  of  n  1  d  considered.  The 
overall  feature  of  R3/R0  in  the  inset  is  dominated  by  the  roughness-scattering  contribution 
Ri/Rq.  The  plateau-like  feature  in  the  range  0.4  x  106  <  Uid  <  0.6  x  106  cm^1  arises  from 
the  abrupt  rise  of  R±  near  n  1  d  =  0.5  x  106  cm _1,  where  ft  passes  through  the  bottom  of  the 
n—  1  sublevel,  where  the  scattering  rate  is  very  large  due  to  the  large  density  of  states. 

The  dominance  of  roughness  scattering  in  Fig.  3  overshadows  the  electron-electron  scat¬ 
tering  as  a  function  of  the  electron  density.  Figure  4  highlights  the  effects  of  electron-electron 
and  phonon  scattering  by  (a)  including  and  (b)  removing  the  roughness  scattering  from  our 
system.  Case  (a)  corresponds  to  the  situation  studied  in  Fig.  3  with  the  same  definitions 
for  Ri,  R2l  f?3  but  not  for  R0.  The  contribution  from  the  phonon  or  electron-electron  scat¬ 
tering  can  be  separated  as  clearly  seen  in  (a)  which  presents  the  scaled  resistance  differences 
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AR/R0  and  in  (b)  which  shows  the  scaled  resistance  R/Ro  as  a  function  of  the  density.  In 
(a),  the  contribution  from  phonon  scattering  is  shown  as  i?2  —  R\  (dashed  curve),  while  the 
contribution  from  electron-electron  scattering  is  shown  as  i?3  —  R2  (solid  curve).  For  phonon 
scattering  in  (a),  i?2  —  R\  displays  a  broad  peak  around  n \ d  =  0.5  x  106  cm-1  where  the 
n  —  1  sublevel  starts  to  be  populated.  Just  above  the  threshold  density  for  the  occupation 
of  the  n  —  1  sublevel,  the  Fermi  wave  vector  kip  of  this  sublevel  is  small  allowing  efficient 
emission  and  absorption  of  small- momentum  phonons.  With  further  increasing  of  nio,  /uf 
becomes  larger,  leading  to  a  suppression  of  phonon  scattering  or  a  decrease  of  R-2  —  R\  with 
n  1  d •  The  electron-electron  scattering  contributions  R3  —  R2  vanishes  in  (a)  at  low  density 
due  to  the  complete  cancellation  of  the  intrasubband  electron- electron  scattering  until  the 
n  —  1  sublevel  is  populated  at  niD  =  0.5  x  106  cm -1.  When  n \ o  increases  above  this  value, 
R/>  —  R2  increases  dramatically  with  uid  because  more  and  more  electrons  are  available  for 
Coulomb  scattering,  in  sharp  contrast  to  the  decrease  of  i?2  —  Ri  with  thd-  As  a  result, 
R3  —  R2  eventually  dominates  R2  —  R\  above  >  0.8  x  106  cm _1.  In  the  absence  of 
roughness  scattering,  we  ford  from  (b)  that  the  effect  of  electron-electron  scattering  (solid 
curve)  is  to  enhance  the  phonon-scattering  resistivity  (dashed  curve)  above  the  threshold 
density  at  n \ n  =  0.5  x  106  cm -1.  A  completely  different  n \ o  dependence  of  the  resistance 
is  seen  in  (b)  for  a  roughness- free  system  in  contrast  with  that  in  the  inset  of  Fig.  3  for  a 
roughness-dominated  system. 


V.  CONCLUSIONS 

In  conclusion,  we  have  investigated  the  effect  of  the  electron-electron  interaction  on  the 
resistance  of  diffusive  electrons  in  a  multi-sublevel  single  quantum  wire  as  a  function  of 
the  temperature  and  the  density  at  low  temperatures.  The  intrasubband  electron-electron 
scattering  is  shown  to  have  no  effect.  However,  for  systems  with  multi- sub  level  occupation, 
electron-electron  scattering  is  shown  to  enhance  the  resistance  significantly.  Also,  this  effect 
is  relatively  more  important  in  high-density  systems.  We  have  demonstrated  that  this  effect 
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can  be  studied  with  arbitrary  accuracy  using  a  formally  exact  solution  of  the  Boltzmann 
transport  equation  in  the  presence  of  roughness  and  phonon  scattering. 
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FIGURES 


T  (K) 


FIG.  1.  Scaled  resistances  {R/Rq)  as  functions  of  the  temperature  T  in  the  range  4  <  T  <8  K. 
Here,  Ro  =  1.06  (h/ 2e2)  is  the  resistance  R\  due  to  surface  roughness  scattering  alone  around 
T  =  0.5  K.  Three  different  cases  are  compared  to  each  other  in  the  figure:  (1)  i?i  (dash-dot-dotted 
curve);  (2)  Roughness  plus  phonon  scattering  (i?2,  dashed  curve);  (3)  Roughness  plus  phonon  plus 
electron-electron  scattering  ( R3 ,  solid  curve).  The  inset  shows  R3/R0  in  the  temperature  range 
0.5  <  T  <  8  K. 
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FIG.  2.  Scaled  resistance  differences  AR/Rq  in  (a)  and  the  scaled  resistance  R/Ro  in  (b)  as 
functions  of  the  temperature  T.  The  roughness  scattering  is  included  in  (a)  but  is  excluded  in  (b). 
Case  (a)  corresponds  to  the  situation  studied  in  Fig.  1  with  the  same  definitions  for  R,\.  R2,  R3 
but  not  for  Rq.  Two  different  cases  are  compared  to  each  other  in  (a)  with  Rq  =  0.106  ( h/2e 2) 
and  L  =  0.1  gm:  (1)  R2  —  Ri  (dashed  curve,  effect  of  phonon  scattering);  (2)  R3  —  R2  (solid  curve, 
effect  of  electron-electron  scattering).  We  have  compared  two  results  in  (b)  for  phonon  scattering 
only  (dashed  curve)  and  for  combined  electron-electron  and  phonon  scattering  (solid  curve)  with 
Rq  =  0.678  (h/2e2)  and  L  =  100  fim. 
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n1D  (106cm') 


FIG.  3.  Scaled  resistances  (R/Ro)  as  a  function  of  the  density  tt-id  in  the  range 
0.8  x  106  <  niD  <  1.0  x  106  cm~1.  Here,  Ro  =  1.236  (h/ 2e2)  is  the  resistance  R\  due  to  sur¬ 
face  roughness  scattering  alone  at  niD  =  106  cm _1.  Three  different  cases  are  compared  to  each 
other  in  the  figure:  (1)  R\  (dash-dot-dotted  curve);  (2)  Roughness  plus  phonon  scattering  (R-2, 
dashed  curve);  (3)  Roughness  plus  phonon  plus  electron-electron  scattering  (R3,  solid  curve).  The 
inset  shows  R3/R0  in  the  range  0.35  x  106  <  nio  <  1.0  x  106  cm'1. 
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FIG.  4.  Scaled  resistance  differences  AR/Ro  in  (a)  and  the  scaled  resistance  R/Ro  in  (b)  as  a 
function  of  the  density  n i d •  Roughness  scattering  is  included  in  (a)  but  is  excluded  in  (b).  Case 
(a)  corresponds  to  the  situation  studied  in  Fig.  3  with  the  same  definitions  for  R,\.  R-2,  R%  but 
not  for  Rq.  Two  different  cases  are  compared  to  each  other  in  (a)  with  Rq  =  0.124  (h/ 2e2)  and 
L  =  0.1  gm:  (1)  i?2  —  Ri  (dashed  curve,  effect  of  phonon  scattering);  (2)  R3  —  R-2  (solid  curve, 
effect  of  electron-electron  scattering).  We  have  compared  two  results  in  (b)  for  phonon  scattering 
only  (dashed  curve)  and  for  combined  electron-electron  and  phonon  scattering  (solid  curve)  with 
Ro  =  0.667  (h/2e2)  and  L  =  100  fim. 
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